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Abstract. We present a detailed calculation of the mech- 
anism by which the Accretion-Ejection Instability can ex- 
tract accretion energy and angular momentum from a 
magnetized disk, and redirect them to its corona. In a 
disk threaded by a poloidal magnetic field of the order 
of equipartition with the gas pressure, the instability is 
composed of a spiral wave (analogous to galactic ones) 
and a Rossby vortex. The mechanism detailed here de- 
scribes how the vortex, twisting the footpoints of field 
lines threading the disk, generates Alfven waves propa- 
gating to the corona. We find that this is a very efficient 
mechanism, providing to the corona (where it could feed 
a jet or a wind) a substantial fraction of the accretion 
energy. 
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1. Introduction 

MHD models have shown that jets can be very efficient 
to carry away the accretion energy and angular momen- 
tum extracted by turbulence from accretion disks (Bland- 
ford and Payne, 1982; Lovelace et at, 1987; Pelletier and 
Pudritz, 1992), if the disk is threaded by a poloidal mag- 
netic field. This fits with the observation that accretion 
and ejection are intimately connected in objects ranging 
from protostellar disks to X-ray binaries and AGNs. How- 
ever these models, based on self-similar analytical com- 
putations or on numerical simulations, most often start at 
the upper surface of the disk. Although more recent works 
(Ferreira and Pelletier, 1993a, 1993b, 1995; Casse and Fer- 
reira, 2000) find solutions connecting continuously the disk 
and the jet, these solutions are heavily constrained by con- 
flicting requirements. These can be traced, in good part, 
to the fact that disk models, whether they rely on specific 
instability mechanisms or on the assumption of a turbu- 
lent viscosity, imply that the accretion energy and angular 
momentum are transported radially outward. They must 
thus somehow be redirected upward to feed the jet. 

The Accretion-Ejection Instability (AEI) of magne- 
tized accretion disks, presented by Tagger and Pellat 
(1999, hereafter TP99), could provide a solution to this 
difficulty. It occurs in the inner region of the disk, in the 
configuration assumed by the MHD models of jets (i.e. a 
disk threaded by a poloidal field of the order of equiparti- 
tion with the gas thermal pressure), and it has the unique 
property that energy and momentum extracted from the 
disk can be emitted vertically as Alfven waves propagating 
along magnetic field lines to the corona of the disk. Thus 
they could provide a source for a jet or a wind formed in 
the corona. 

This ability to emit the energy and angular momentum 
as Alfven wavewas recognized in TP99, and indeed justi- 
fied the name given to the instability. However this pos- 
sibility was shown only in a WKB approximation, valid 
away from the region (the corotation radius, where the 
wave rotates at the same velocity as the gas) where Alfven 
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wave emission is most efficient. The WKB result was found 
divergent at corotation, providing a good indication that 
this mechanism of vertical emission could be quite effi- 
cient. 

The goal of this paper is to present a more general 
derivation, valid in the corotation region. We use a de- 
scription of the waves in three dimensions (whereas TP99 
was basically a model averaged over the disk thickness). 
We are thus able to give an explicit computation of the 
Alfven wave emission mechanism and of its efficiency. The 
main unknown to be solved for at this stage is the fraction 
of the accretion energy and angular momentum, extracted 
from the disk by the instability, which will end up emitted 
to the corona. 

The result we present is quite limited: the constraints 
of giving an analytical derivation force us to use a very 
simplified magnetic field geometry, namely an initially 
constant and vertical field, and a more realistic one would 
certainly affect the result. Appendix [A| is dedicated to 
the case of a radially varying B z field. We show that the 
present computation may be apply in the case of a slowly 
variable B field. On the other hand, the result we obtain 
is interesting in itself: we will show that, in linear theory, 
the flux of the Alfven waves is again divergent at corota- 
tion. Although we discuss how it can be regularized, our 
interpretation will be that the efficiency of the mechanism 
is indeed quite high, but that we are reaching the limits 
of linear theory and that the true result will most cer- 
tainly be determined by self-consistent non-linear effects. 
The present work should thus be viewed as an exploration 
of the complex physics involved and of its potential effi- 
ciency, which will then have to be treated by non-linear 
simulations. The paper is structured as follows: in the next 
section we will briefly review the main properties of the 
instability, and its interest to explain the low-frequency 
Quasi-Periodic Oscillation (QPO) of X-ray binaries. Sec- 
tion U will present the system of equations to be solved, 
and section ^ their combination into a variational form 
containing the physics of the problem. In section [5] we will 
compute the Alfven wave flux, and we will discuss the sig- 
nificance of this result in section ^. In appendix ^ we will 
present the variational form we obtain in a more general 



Varniere and Tagger: AEI: Feeding the Corona with Alfven Waves 



3 




Fig. 1. The structure of the instability is shown here schematically as a function of radius. It is formed of a standing 
spiral density wave in the inner part of the disk, coupled to a Rossby vortex it excites at its corotation radius. The 
spiral grows by extracting energy and angular momentum from the disk, and depositing them in the Rossby vortex; 
the latter in turn generates Alfven waves propagating toward the corona of the disk. 



geometry with a radially varying B z and the limitation it 
implies. 



2. The Accretion-Ejection Instability 

We will give here a short review of the main characteristics 
of the AEI, and refer the interested reader to TP99, or to 
Varniere et al. (2002) and Rodriguez et al. (2002) where a 
detailed comparison with the low-frequency QPO of black- 
hole binaries is given. Non-linear MHD simulations were 
performed by Caunt and Tagger (2001). The AEI is es- 
sentially a spiral instability similar to galactic ones but 
driven by magnetic stresses rather than self-gravity. It af- 
fects the inner region of an accretion disk threaded by a 
poloidal magnetic field of the order of cquipartition with 
the thermal pressure of the gas {(3 = 8np/B 2 of the or- 
der of unity), i.e. the configuration and magnitude used 
in most MHD jet models. 



The instability is composed of a spiral density wave 
and a Rossby wave it generates at its corotation radius 
(the radius where the wave rotates at the same velocity as 
the gas) . The spiral forms a standing pattern between the 
inner edge of the disk and the corotation radius. Because 
of differential rotation, it couples to the Rossby wave, i.e. 
the spiral and Rossby waves exchange energy and angu- 
lar momentum in the corotation region where, because of 
differential rotation they loose their separate identities. 

In this process the Rossby wave also forms a standing 
vortex, and the spiral grows by storing in it the energy 
and angular momentum it extracts from the inner region 
of the disk (thus causing accretion). 

In a thin disk in vacuum, the process stops there. This 
is the case in galaxies and recently Fridman et al. (2002) 
observe and analyse such a vortex in NGC 157. However 
one must remember that the Rossby vortex represents a 
torsion of the footpoints of the field lines threading the 
disk. Thus if the disk is covered by a low-density corona, 



4 



Varniere and Tagger: AEI: Feeding the Corona with Alfven Waves 



this torsion will propagate vertically along the field lines as 
Alfven waves: these will thus in turn transport a fraction 
of the accretion energy and angular momentum to the 
corona, where they could provide a source for a jet or 
an outflow. The whole process is shown schematically in 
figure |l|. 

3. Linearized equations 

3.1. Unperturbed equilibrium 

As in previous works, we consider the very simple 
setup of an infinitely thin disk threaded by a vertical 
magnetic field adding this time the hypothesis that it 
is radially constant. Here however the disk is embedded 
in a low density corona. This simplified geometry will 
be enough to fully characterize the emission of Alfven 
waves in the corona. The equilibrium field, B° = Boe z , 
is assumed to be of the order of equipartition with the 
gas pressure (plasma (3 = 8itp/Bq ~ 1) in the disk. The 
case of a radially varying Bq will be studied in appendix [A]. 

We choose to present here results with a constant B° 
because assuming that B° depends on r creates at equi- 
librium a radial magnetic pressure force which can easily 
become dominant in the corona, since the other forces 
(gravity, Coriolis, pressure) act on a very low density. As 
a result a realistic equilibrium will in general include a 
flow along the field lines, and lead to a complex configura- 
tion limiting our ability to extract analytical results unless 
artificial assumptions are made. 

On the other hand instability requires, as found in 
TP99: 

where W = n 2 /2n and n — fl ~ r 3 / 2 in a Keplerian disk. 
Thus, for this quantity to grow with r, our equilibrium 
model needs to have a disk surface density growing with r, 
fast enough for this condition to be fulfilled. This appears 
as an ad hoc model which does not claim to represent 
a real diskjj]. We emphasize however that it is physically 
consistent, so that our computation of the Alfven wave 
flux remains exact and will allow us to explore physics 
which applies to more general disk models. 

3.2. Pertubations 

We work in cylindrical coordinates [(s = lnr),6,z\. We 
consider linear perturbations, described by Lagrangian 
displacements £ with: 



1 It is worth noticing (J. Ferreira, private communication) 
that the self-similar assumption used in MHD models of jets, 
following Blandford and Payne, always results in this deriva- 
tive being positive and equal to +1/2, as a result of the self- 
similarity laws. Thus all these models are unstable to the AEI. 



where V is the perturbed velocity, Vq = rfl(r) is the 
equilibrium rotation velocity, and 
d d 

In TP99 the instability was studied by giving a full solu- 
tion for its spatial structure in cylindrical geometry; the 
emission of Alfven waves was then computed perturba- 
tively, far from the corotation radius, in the WKB limit. 
The waves were described by their compressional and tor- 
sional components, 

D = Vj_ ■ V± and R = V_l x V± 

where the subscript _L describes the components of a vec- 
tor in the plane of the disk. We find this description in- 
convenient here and rather write the displacement as: 
Bo(r)£± = -V_l$ - e z x 

where again $ and 4" represent the compressional and 

torsional parts of the displacement. 

From this we write the induction equation, 

d t B + V x(V xB) = 0: 

giving the perturbed magnetic field as: 

Bi = (B°-V)£ x =B cUx (1) 
and from V ■ B\ = 0: 

Bl = -V±-(B Z ± ) (2) 
Comparing with the solution of the continuity equation, 
which gives the perturbed density p = — Vj_ • (po one 
notes that equation || is an equation of conservation for 
the vertical magnetic flux. 

We consider perturbations varying as exp[i{md — cut)] 
(so that m will be the number of arms of the spiral) . We 
define a de-dimensionalized perpendicular gradient, 
V±X = rViI = d s Xe r + imXe$. 
We finally get the perturbed magnetic field: 
imd z ty 



B 1 = 



I 



-imd z $ - d ZjS ^ 



(3) 



r -L 

As in TP99, our goal is to derive from the MHD 
equations a variational form. Its real (self-adjoint) part 
represents the energy of the perturbations and describes 
their structure in the disk. Amplification and damping 
appear through imaginary terms, representing a flux 
of energy: either at the boundaries of the system (as 
outgoing waves) or at the corotation resonance. The 
classical physics of spiral waves in disks with differential 
rotation is that they have positive energy beyond the 
corotation resonance (rc such that m£l(rc) = to), and 
negative energy inside it: at r < rc the gas rotates faster 
than the wave, whose presence decreases the total energy 
of the system (it releases more gravitational energy than 
it costs kinetic energy). 

The spiral can thus be amplified if positive energy is 
emitted beyond corotation (as a spiral wave emitted 
outward: this is the Swing mechanism, responsible for the 
amplification of galactic spirals), or stored in a Rossby 
vortex at the corotation radius: in TP99 this was shown 
to result, in a magnetized disk, in the Accretion-Ejection 
Instability. This corotation resonance introduces, in 
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the variational form of TP99, a pole (a denominator 
proportional to [lj — mO(r)], singular at corotation). 
This pole contributes, following the Landau prescription 
familiar in plasma physics, an imaginary term which 
represents the energy exchanged with the Rossby vortex. 



concerns the slow wave, which should thus not be so 
important for us. We thus postpone its consideration 
to future work. In the same manner, consistent with 
the thin disk approximation, wc will neglect the vertical 
component of the Euler equation. 



In TP99 it was also found that, when one takes into 
account a low-density corona above the disk, Alfven 
waves are emitted toward the corona; their contribution 
was computed in a WKB approximation, valid only away 
from corotation. This had of course a limited interest, 
since we expect the Alfven waves to be strongest at 
corotation where the torsional motion associated with the 
vortex is strongest. Here we will present a more general 
formulation, showing how the energy flux of the Alfven 
waves appears as an additional imaginary term in the 
variational form, generated at the vertical boundaries 
of the disk. This will allow us to compute explicitly the 
Alfven wave flux at corotation. 

In terms of the Lagrangian displacement, the linearized 

can 



i(m$— u)t) 



Euler equation for a perturbation varying as e 
be written as: 

(-w 2 + 200') i r + 2iCjQ f„ = f r (4) 
— lu 2 ^ — 2iojQ £ r = (5) 

-u% = /, (6) 
where / = F/po, F is the force acting on the system, 
po(s,z) is the equilibrium density, uj(s) = lu — mfi(s), 
and the prime denotes differentiation with respect to s. 
In terms of $ and 4" the r and i? components of these 
equations become: 



-.2 



\0J 



200') d s $ + 2mcDO$ 



-2iQild„V - im{Co 2 - 200') * = rB f r (7) 
2«a)09 s $ + «mw 2 $ + w 2 a s * + 2ma)0* = rB fo (8) 
For simplicity we will neglect the pressure stress, since 
it was found in TP99 to be negligible, compared to the 
magnetic ones, in the physics of the instability. It could 
be added without change on the emission of Alfven 
waves. On the other hand, this means that we do not 
consider here the slow magnetosonic wave which could 
also propagate upward from the disk. We have two 
reasons for that: the first one is that this would lead to 
a very complex problem, which goes beyond the goals 
of the present paper; it would force us to consider the 
detailed physics at the slow magnetosonic point above 
the disk surface (see e.g. Ferreira and Pelletier, 1995, 
and references therein), where the gas is first extracted 
from the disk. This means that, in stationary MHD 
models of jets, the physics of the slow magnetosonic 
wave concerns the mass loading of the field lines, rather 
than the acceleration of the jet which is more associated 
with the Alfven wave. Our second reason is that we will 
show below how the Alfven wave is mainly excited in a 
resonant layer at the corotation radius, by the singular 
Rossby vortex generated in the disk. No such singularity 



Linearizing the contribution of magnetic stresses, (j x 
B)/p, we find : 
j 1 x B° 



f = 



Po 



(9) 



where the subscripts and 1 note equilibrium and per- 
turbed values. We thus get the expression of the magnetic 
stresses: 

B fr = Vl{d z Bl - d r Bl) 

BoU = Vl{d z Bl - d 6 B\) 

(where V\(z) — Bq/Attpq); we express this in terms of $ 
and 'F, and get from equations (@-|(): 
rB Q [(-Co 2 + 200') £ r + 2icD0 = 

-Vl{d s V 2 <S> - imd z 2^) (10) 
rB [-Cj 2 & - 2iQQ £ r ] = -Vl{d St ^ + imV 2 f) (11) 



4. The system 

In order to obtain the new variational form, we first write 
a quadratic form built from the divergence and the curl 
of \Bq times the Euler equation, leaving A unspecified for 
the moment. 

We first apply the operator r 2 Vj_ • XBq to the Euler equa- 
tions, i.e. we compute: 

^(rAB @) +imr\B ^ 
we get: 

V ± (A£ 2 V ± $) - 2a j (AOO'<9 s $) 
+ 2to9 s (AOw)$ - 2iV ± (AwOV_L*) 
- imd s {\Co 2 )^ + 2imd s (\Q.Q,'^) 
= -AV^ViV 2 $ 

-{XdsVl + Vp s \)(imd z 2V - 9 S V 2 $) (12) 
Then applying the operator r 2 V^ x \B§ to the Euler 
equations: 

V±(A£ 2 V±*) + 2rnw<9 s (AO)>F 
+ 2iVj_(A£0V_L$) + imd s (XQ 2 )^ 
+ 2imAOO'<9 s $ 

-(Xd s Vl + Vld s X){imd z 2^ - 9 S V 2 $) (13) 
In order to get a variational form from these equations we 
need to get rid of the term in (imd z 2^— <9 S V 2 $), present in 
both. This will be the case if A9 s Vf + V|9 S A = 0. We thus 
make the obvious choice A = V^ 2 (s, z), and equations ( |l^ - 
[ilt) then become: 

'00' 



V 2 A 
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, flu) , 
+2md s | ^2 | <I> - 2/Vj 



-imd s [ ^ I * + 2im3 s 



V. 



n 



and: 

A 



(14) 



= -a z2 Vi* (15) 
We will later put these equations in a more compact form, 
but use for the moment the present one which is best 
adapted to a variational formulation. 

4-1- The Variational Form 

In order to obtain the variational form we will first inte- 
grate these equations vertically. Assuming that the disk 
is covered by a low-density corona, we integrate between 
z = — z maz and +z max , chosen well into the corona so 
that at that height only the Alfven 2 waves propagate and 
their vertical propagation can be described in the WKB 
approximation. We write the integral: 



F 



$* [14] + * p|* dzds 



= (16) 
where |ll|* is the complex conjugate of equation (|l5|), 
and the radial boundaries s m i n and s max will be discussed 
later. We will find that, thanks to our choice of A, this 
form has the properties we wished: it is composed of 
a main part, which is variational, and additional terms 
which can be treated perturbatively to give amplification 
or damping of the waves. 

After some algebra, integrating by parts and grouping 
terms we get: 



-S (|V ± *| 2 -|Vx*| 2 ) 

V A 

n 

2m d s (^) |$| 2 
2mu>d.(J^\ *| 2 



$*Vi V 2 $ Us dz 



2 As shown in Tagger et al. 1990 the fast magnetosonic wave 
does not propagate in the corotation region. 



ax dz 



Co 1 $*V_l$ - 2fMT $*<9 S $ 



w 2 *Vj_** + 2i m nn' 

2iQfl ( <JV_l$* -I- $*V_l* 



ds 



(17) 



Equation (17) is equivalent to the variational form de- 
rived in TP99, but the vertical integration (rather than the 
approximation of an infinitely thin disk) will give us ac- 
cess to the emission of Alfven waves. The first five terms of 
equation^]) appear to be obviously hermitian, although 
they hide the imaginary contribution from the corotation 
resonance (i.e. the growth or damping term corresponding 
to the energy exchanged with the Rossby vortex) which 
will be discussed in the next subsection. 
The right-hand side is formed of boundary terms, ob- 
tained in the integrations by parts. All these terms are 
easily shown to be imaginary (i.e. correspond to growth 
or damping) , if the boundaries are far enough that the ra- 
dial and vertical derivatives can be estimated in a WKB 
approximation, d s = ik s , d z = ik z , and if waves do prop- 
agate at the boundaries, i.e. if the wavevectors are real so 
that waves can effectively transport energy away. 
The first bracket corresponds to the flux of the wave at the 
radial boundaries; as explained in TP99, a wide range of 
boundary conditions at the inner disk edge allow the waves 
to be reflected without loss of energy, i.e. this term van- 
ishes or remains real at s m i n . At s max it gives the flux of 
an outgoing wave, responsible for the usual Swing amplifi- 
cation of spiral waves (driven by self-gravity in galaxies, by 
pressure in the Papaloizou-Pringle instability (Papaloizou 
and Pringle, 1985), or by magnetic stresses in Tagger et 
al. 1990). 

The last term is new and corresponds to a flux at the ver- 
tical boundaries, i.e. the flux of the Alfven waves emitted 
vertically. This is confirmed by the fact that this term is 
associated with the torsional ('J) component of the per- 
turbations. 



4-2. Corotation Resonance 

In our variational form, equation (|l7|), the corotation res- 
onance does not appear explicitly as it did in the equiva- 
lent form of TP99: we do not get denominators containing 
u>, vanishing at corotation. But the corotation resonance 
must of course be present, since the physics described here 
is more general than in TP99: it is hidden here in the sin- 
gularity of \I/: equation (Jl^) has a singular point at coro- 
tation (where ui vanishes) because the highest-order radial 
derivative of \& is proportional to w 2 (assuming, as will be 
checked a posteriori, that the vertical derivative vanishes 
at corotation), while the terms in $ are proportional to uj. 
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We thus turn, in the vicinity of corotation, to a Frobenius 
expansion of the form: 

* = a_iw _1 + a n uj n + lncD + b n Jj n baa) (18) 

$ = c + c„a)" + d\Cb In w + tf„a>" In a) (19) 
where the coefficients a, b, c, d depend on z. We first 
note that, in equations ( p^]l5| ), the force terms (in the 
right-hand sides) are small in the disk, since they are 
of order k 2 V^/Vl 2 (where k is a radial or vertical scale 
length) compared to the other ones. Since we consider 
that (3 — 2c 2 s jV\ ~ 1, and c s ~ hCl/r where h is the disk 
thickness, these terms are small near corotation (they play 
of course their role in the radial propagation of the waves, 
far from corotation, see TP99). The vertical wavenumbers 
will be discussed in more details in section |5.l| below. We 
re-arrange equations (|l^-|l5|) as: 



Vi 



Co 2 - 2nn' - ^ 



2mJ.iK I I <T> - 4/// 



4> 



-2iV 



v, 



2 Vl(u$) 



and: 



imd s | ) * 



n 



2n>l-(K I -pr I 



Q - 



n 



V i $ 



= -a 22 vi* (21) 

and neglect in the disk the force terms in the right-hand 
sides. To leading order (a) -1 ) equation ( ^p|) gives: 

di = ~b (22) 
while to order Jj° equation (|2l]) gives: 

w\ n' 2 , nn' 2 



so that 



ttt + 2imd\ 







^rrC'sln 



w 



(23) 



(24) 



and the following coefficients can be derived from the next 
orders. 

The manner in which $ and '5 project on the solutions 
which are regular and singular at corotation depends on 
the global solution, which must be obtained numerically 
as in TP99, of the problem with its boundary conditions. 
Here we will only use the fact that \& has a singular con- 
tribution (whose Frobenius expansion starts with a_io) _1 ) 
which makes the quadratic form, equation (^) , non vari- 
ational at corotation. Combining the singular terms, we 
find their contribution to the integral: 



Ft 



Corot. 



Jds Jdz [ ^ |d s <F 



2mu>d fS I — w 



(25) 



Integrating by parts the first term, and retaining only the 



term in the expansion of 'F, we find: 



Corot. 



ds 



dz M-a 



ds — o 



\M> (26) 

where S is the surface density of the disk. As expected we 
recover the result from TP99 that the resonant contribu- 
tion to the variational form (and hence to the growth rate 
of the instability) is proportional to d s {pW/ Bq), and will 
give an imaginary contribution from the pole at corota- 
tion. 

Our goal now is to compute how this singular vorticity 
generates Alfven waves, transmitting to the corona a part 
of the accretion energy and angular momentum extracted 
from the disk. This flux is readily identified, in equation 
( |l7| ) , as the last term which represents a contribution from 
the lower and upper boundaries of our integration domain: 
this is thus the flux emitted vertically, and we expect that 
the large amplitude and the singularity of 'F will give a 
strong contribution at corotation. 



(20) 5. Dispersion relation and Alfven Flux 

5.1. Dispersion Relation in the Corona 

The flux of Alfven waves to the corona appears in our 
variational form, equation ( |l7| ) , as a surface term taken at 
the lower and upper boundaries of our integration domain: 



F Alfven 



<F* d z Vi * 



(27) 



We take these boundaries far enough above the disk that 
the density varies smoothly, so that the vertical derivative 
can be computed in a WKB approximation, d z ~ ik z . A 
WKB approximation in the radial direction was also used 
in TP99 to compute this flux. However, near corotation, 
the singularity of >F makes this radial WKB expansion 
invalid. We will show this by turning again to a Frobenius 
expansion, but retaining the force terms in equations ( pOf - 

Let us first consider the case of a disk in vacuum: 
the Alfven velocity in the corona goes to infinity, so that 
equations ( 20 2l]) reduce to A$ = and <9 Z 2>F = 0. In 
a radial WKB approximation, this gives the result famil- 
iar in spiral wave theory, that $ varies above the disk as 
exp(— \kz\), where k is the horizontal wavenumber, so that 
$ vanishes exponentially. On the other hand one finds that 
^ stays constant with z. The full vertical solution, valid 
across the disk, was given by Tagger et ah, 1992. If the 
coronal density is now small but non vanishing, equation 
( pp| ) shows that far above the disk $ must be of the order 
of (r2 2 /(fc 2 V^4 00 )^>, where fc _1 is the radial scalelength of 
$, and Vaoo is the Alfven velocity in the corona, large but 
not infinite. Thus now $ is negligible in equation (pl|), 
while the vertical derivative of 'F must be retained. In a 
WKB approximation in z we now expand this derivative 
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as: 

d z 2 « -k\ = -e 2 u) 2 - e 3 ^ 3 + . . . (28) 
and to lowest order of the Frobenius expansion equation 
( ^l| ) now gives e 2 — and: 

/ W \ si' 2 

a -i d s I 775— I + TO& o— 2e 3 mfi' 2 a_i=0 (29) 



v Ac 



v Ac 



Since the result obtained at the disk, equation (p4[), deter- 
mines the values of a_i and & & t the base of the corona, 
we can use them to get giving finally (assuming for 
simplicity that the radial profile of V\ is independent of 

z): 

ktvL 



QQ S (W 



(30) 



This result comes as a surprise, for two reasons: the 
first one is that, since we consider the propagation of 
a perturbation which is very localized in s (since it 
is singular), one would have been tempted to use a 
radial WKB approximation to solve equation (|2l]). From 
the first terms in the left and right-hand sides of this 
equation, one would get k' l V\ 00 = a) 2 , i.e. the familiar 
dispersion relation of Alfven waves, Doppler-shifted by 
differential rotation. The answer here is simple, since 



with ^ 



the first term in the right-hand side of 



equation (|2lj) vanishes to lowest order in oj~ l . Thus the 
'usual' dispersion relation applies to regular features, but 
not to the singular one we are concerned with. 

The second surprise is that k 2 is proportional to d> 3 , 
so that we find k z real only on one side of corotation, de- 
pending on the sign of the derivative of WjV\. Thus the 
singular perturbation will propagate only on one side of 
corotation! This can be understood by returning to equa- 
tion and writing it now in terms of Y = (D\l/, still 
neglecting the contribution of <&. We get: 

u>V±^S*±Y + mYd s = klV\ (Jpj (31) 

Let us first neglect the magnetic tension force, on the 
right-hand side. Equation ( plf ) gives, in a radial WKB 
approximation (d s w ik): 
W\ 

(32) 



raV\ 



n 



k 2 + m 2 

This is just the dispersion relation of Rossby waves, 
Doppler-shifted by differential rotation and including the 
radial gradient of V%\ Because this dispersion relation is 
odd in ui (whereas more usual ones are even) it allows 
Rossby waves to propagate only on one side of corotation^. 



The meaning of equation (31) becomes simple: it de- 
scribes how Rossby waves, forming the singular part of 
our full solution for 'J', will now propagate upward along 



For a more detailed discussion of Rossby waves with differ- 
ential rotation, see Tagger, 2001. In particular, this shows how 
a Rossby wave in a disk always collapses to a singularity, the 
one analyzed here, after a finite time, of the order of a rotation 
period. 



the field lines as what we might call Rossby- Alfven waves. 
If the radial derivative of WjV\ is positive (this is the 
instability criterion of TP99), the wave propagates only 
beyond corotation and thus carries a positive energy flux, 
so that its formation in the disk and propagation along 
the field lines destabilize the negative energy spiral wave 
inside corotation, which is the main component of our in- 
stability. 

This makes the physics of the coupling between the insta- 
bility and Alfven waves much more complex than the res- 
onance, found by Curry and Pudritz (1996), between nor- 
mal modes of the Magneto-Rotational Instability (MRI) 
and Alfven waves. The main reason is that they work in 
a vertical WKB approximation, assuming that the modes 
have a fixed vertical wavenumber k z ; it is worth remem- 
bering here the result of a more complete vertical solution 
(Tagger et ai, 1992): the AEI corresponds to solutions 
with n z = nodes across the disk height, whereas the MRI 
corresponds to solutions with n z < 1. For n z 3> 1, most 
unstable when 3> 1, a vertical WKB approximation can 
be used. Keeping k z fixed results in finding a resonance 
away from corotation, at Q = k z VA, and thus ignoring 
the coupling with the Rossby vortex. It also ignores the 
difficulties, encountered here, associated with the vertical 
variation of the Alfven velocity 

5.2. Alfven waves Flux 

We can now compute the flux of Alfven waves, appearing 



in the right-hand side of equation (17) 



F Alfven — i 



[MVx*| 



ds 



(33) 



where we have assumed that the upper boundary is high 
enough above the disk, so that the vertical derivative can 
be replaced by ik z . This contribution was computed in 
TP99, in a radial WKB approximation valid far from 
corotation. Here our goal is to compute the localized con- 
tribution of the resonance, i.e. the flux of Alfven waves 
launched from the Rossby vortex, and more precisely the 
flux associated with the singular part of \& as computed 
in the previous sections. Thus we retain ^ rj a_iw _1 near 
corotation, and take into account only the region beyond 
corotation where (if the radial gradient of W jV\ is pos- 
itive) the Rossby-Alfven wave propagates, so that k z is 
real. We choose for k z the sign of z, in order to have waves 
propagating away at z min and z max . We get: 



F 



Alfven 



2i\a. 



ds 



LU 3 / 2 



v A 



w 



"I 1/2 



(34) 



This integral is strongly divergent at corotation, as ex- 
pected. Various effects can regularize it, in particular the 
presence of an imaginary part of u>, or the fact that most 
of our derivations fail at radial scales of the order of the 
disk half-thickness h. Let us consider in particular equa- 
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tion (|20|): we have solved it in the disk by a Frobenius 
expansion, neglecting altogether the force term, on the 
right-hand side, although it contains the highest-order ra- 
dial derivative. Our reason was that this term is small as 
long as Co is not too small: remembering that we assume 
the plasma (3 to be of the order of 1 in the disk, we have: 



Vi 



4 



r 



(35) 



so that this term, acting on the divergent part of is of 
the order of 



n 



r j \w ; 

compared to the one we retain. This means that our ex- 
pansion breaks down when uj ~ (h/r) Q, or s ~ h/r, and 
we will take this as the lower bound in our integral, equa- 
tion (|34|). It is quite possible that the full solution, at 
s < h/r, would still give a divergent result, but then the 
growth rate of the instability, estimated in TP99 to be 

Im(v) ~ (^) Re(uj) 

would give a similar regularization of the integral. Thus 
we get a conservative estimate for the flux transported 
toward the corona by Alfven waves: assuming that all the 
gradients are of the order of 1, we get: 

F Mfven ^\a^ ^ (36) 

On the other hand, the flux deposited at corotation in 
the vortex (which is the energy removed from the central 
region of the disk, causing accretion) is given by equation 
(p6|). As in the classical problem of Landau damping, its 
imaginary part is given by a Cauchy residue: 



Corot. 



a 



{Y,W\ 



(37) 



m\fl'\ ~°\Bl J 
Both fluxes are related to the singularity of '5 at corota- 
tion, and thus proportional to |a_i| 2 . Assuming again for 
simplicity that (3 ~ 1 and that the radial gradients con- 
tribute numbers of order unity, we thus get for their ratio 
a very simple expression: 



(38) 



F Alfven ( Poo\ /^\ 3/2 
Fcorot. \PD ) \h) 

where we have used £ ~ prjes/Q, and po, Poo are the 
densities in the disk and in the corona. The first term 
in the right-hand side is a priori small. It is typical of 
magnetic braking processes, and leads to a weak efficiency 
of mechanisms of magnetic origin coupling the disk to 
the corona. The good surprise for us here is that it is 
multiplied by the second term, which can be quite large 
(typically a few hundred to a thousand in the disks of 
X-ray binaries). The flux ratio can then be a significant 
fraction of unity, as soon as the corona has a density 
which is not vanishingly small. For higher density making 
our estimate larger than 1, a full 3D computation taking 
into account the finite thickness of the disk would be 
required. 



5.3. Non-Vertical Unperturbed Magnetic Field 

In the body of the paper we have studied the case of a 
vertical and constant unperturbed field. In Appendix ^ 
we discuss the effect of a field which is still straight but 
depends on the radius. The analytical computation cannot 
be fully completed, mainly because no simple equilibrium 
exists without a flow along the field lines; but we show that 
within reasonable bounds (that the magnetic stress term 
does not exceed the centrifugal and gravitation forces) our 
conclusions should not change. 

On the other hand a realistic model should also include 
the curved geometry of the magnetic field, as obtained in 
jet models. In this case the situation becomes much more 
complex because in general the three basic MHD waves be- 
come coupled by the geometry, in addition to the coupling 
by differential rotation studied here. In a straight field the 
Alfven and slow magnetosonic wave are decoupled: this 
has allowed us to defer the consideration of slow magne- 
tosonic waves, together with vertical motion, to separate 
work. In a curved field the Alfven wave will include ver- 
tical motion, and we should in principle include all three 
components of the displacement and all three MHD waves. 
The mixing can still be weak if the scales are very differ- 
ent. For instance, the characteristic wavelength of the slow 
wave is of the order of cg/iD, i.e. of the disk scale height 
if we are not too close to corotation. Thus if the mag- 
netic field is curved on a large scale (of the order of r), 
the mixing of the waves is weak and our conclusion should 
not change much. In particular the Alfven wave is excited 
much more efficiently than the slow wave, for which we 
expect no singular source analogous to the Rossby vortex 
for the Alfven wave. 

On the other hand more elaborate jet models (e.g. 
Casse & Ferreira, 2000) show that a slow magnetosonic 
point forms above the disk. The field lines are sharply 
bent in its vicinity. But there the Alfven velocity is al- 
ready much higher than the sound velocity. This dispar- 
ity should again maintain a separation between slow and 
Alfven waves, and allow the latter to propagate the vortic- 
ity from the Rossby wave. A realistic computation goes far 
beyond our present abilities. But we note that the coupling 
of waves by geometric effects might introduce interesting 
new channels to deposit energy and momentum from the 
wave in the corona. 



6. Conclusion 

We have presented a computation of the flux of Alfven 
waves emitted to the corona of a magnetized disk by the 
Accretion-Ejection Instability. This means that we have 
justified here the name chosen by TP99: the instability 
is a spiral density wave, which grows by extracting en- 
ergy and angular momentum from the disk (thus causing 
accretion) and transferring them radially outward to the 
Rossby vortex at corotation; a significant fraction, given 
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by equation (p8|), of this flux is then transmitted verti- 
cally to the corona as Alfven waves. We expect that, if 
the Alfven waves can deposit their energy and momen- 
tum in the corona, this would be an ideal mechanism to 
feed a wind or jet directly from the accretion process in 
the disk. 

The amplification of the wave (and thus the flux de- 
posited by the spiral in the vortex) and the flux transmit- 
ted to Alfven waves are both linked to the singularity of 
the vortex. This allows us to give in a very simple form a 
result of paramount importance in the physics of accretion 
disks and jets: an estimate of the fraction of the accretion 
energy, extracted from the inner region of the disk, which 
will end up in the corona where it might feed a jet. This 
fraction is of the order of unity if the coronal density is not 
too low (typically 10~ 4 of the density in the disk would 
be sufficient, in an X-ray binary). 

In order to obtain analytically a physically consistent 
result, we have had to use a very artificial configuration of 
the equilibrium magnetic field, vertical and independent 
of r. On the other hand this has allowed us, proceeding 
rigorously by perturbation of a variational form, to obtain 
an exact result clarifying the role and the physical nature 
of the singularity of the Rossby vortex at corotation. We 
can thus expect that these results would survive less strin- 
gent assumptions on the equilibrium configuration. 

However this must be taken carefully: our final result is 
in fact divergent at the corotation radius, and regularized 
by the effect of the finite thickness of the disk, or by the 
growth rate of the instability. In both cases, it depends on 
the density in the corona of the disk. Thus we believe that 
the end result should be obtained from a self-consistent, 
non-linear description where the growth of the instability 
itself affects the evolution of the magnetic geometry and 
the mass loading of the corona. 

In this respect it is worth mentioning one of the results 
of stationary MHD jets models: in these models, as the gas 
is accelerated along the field lines it passes a slow magne- 
tosonic point where the field lines bend outwards. Magne- 
tocentrifugal acceleration can then proceed and leads to 
the formation, higher up and further out, of an alfvenic 
point. The slow magnetosonic point is thus associated with 
the mass loading of the field lines, and the alfvenic point 
to the acceleration. By analogy we can thus expect that, 
while the Alfven waves described in the present work al- 
low to accelerate the gas, the instability can also generate 
slow magnetosonic waves which will lift the gas above the 
disk. The coupling of the instability to the slow wave will 
be the object of a forthcoming paper. 

Appendix A: Radial dependence of B z 

In the body of the present article we made the simple 
assumption that the equilibrium field B° was vertical 
and constant, allowing us to get an analytical derivation 
of the Alfven waves flux. In this appendix we give the 



general derivation in the case where the field is still 
straight along z but depends on s. We will follow the 
same computation as in the main text, referring to the 
corresponding equations. 

The first modification appears in the contribution of 
magnetic stresses (j x B)/p. We now have to take into 
account the equilibrium current, j : 



f = 



j" x B 1 



x B° 



- if x B°) * 
Po 



(A.l) 



Po Po 
We get pi from the continuity equation, 

Pi = -V • (po€) (A.2) 
which gives after some transformations, and using equa- 
tion (|): 

Bo— = Bl-\g{rB ^ r ) (A.3) 

Po r 2 
Where the function g is defined by: 

d 2 



9(rB £ r ) 



Po_ 

Bn 



0J 2 p 



[Bo] 



0:. 



(rB ti r ) 



(A.4) 



and rBa£ r — — d s 3> + imty . Here angular brackets note the 
logarithmic derivative with respect to s, [X] = 1/X d s X, 
and g comes from the departure, due to magnetic pressure, 
from the keplerian rotation curve. It is weak in the disk 
for /3 ~ 1; it might become important in the corona, but 
this would require to take into account the velocity flow 
associated with the formation of the jet. This full compu- 
tation is beyond our present capabilities. 
We thus get the expression of the magnetic stresses: 

d r Bl) — B Q — — -Bl + B — — -B — 

Po Po Po 

Boh = Vl{d z Bl - dvBl) 

Expressing this in terms of $ and ^ we get: 

rB [{-Co 2 + 2nn') £r + 2iu>Q &>] = 



Bofr = Vl{d z Bl 



-Vl{d s V 2 § - imd z 2^) 



Bq [— Co 2 ^ — 2iuj£l £ r ] 



^ [B Q }g(rB ^) 



(A.5) 



-Vl{d s ^ +imV 2 $) (A.6) 
which generalize equations and 

A.l. System 

From these new equations, using the same method as in 
section ^ we obtain the following parametric system: 
V ± (Aw 2 Vj_$) - 2d s (\nn'd s <i>) 

+ 2md s (AQ£D)$-2iVx(Aa>fiV x #) 
- imdsiXu 2 )^ + 2imd s {\Q.Q!^) 

= -xvlvlv 2 ® 

-{Xd s Vl + Vld s X){imB 7 2^> - 9 S V 2 $) 



-(av; 



1 1 \d s V\ 



Vp s \)d s { [ ^Lg{rB^ r ) ) (A.7) 



and 



V ± (A£ 2 V ± *) + 2mQd s {\Q)y 
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+ 2iV _l(AwOV + imd s (\Cj 2 )§ 
+ 2im\nn'd s <f> 
= -AV^Vi* 

-{\d s Vl + V2d s X){imd z 2^ - 9 5 V 2 $) 

+ (AVj + \d s Vl + Vp s \) l -^[B Q } g{rB^ r ) (A.8) 

In order to get rid of the term in (imd z 2^ — d s V 2 $) we 
make the same choice as in the constant- B° case, namely 
A = V^ 2 {s, z). The system equivalent to ( p^]l5| ) is: 



V ( —,V <I> ] - 2,7 



, w 2 \ /fifi' 
I ^2 J * + 2imds y~yT^ 



(A.9) 



and: 



+2jm-r7^-9 s $ 



A. 2. The Variational Form 



(A.10) 



After some algebra, integrating by parts and grouping 
terms we get the equivalent of equation (fl7|) 



$ (|Vl*| 2 

V A 
"A 

2m & ['^] \9\ 



|V,4>| 



2mw 9 S 



f2 

$*V 2 i v 2 $ 



1*1 



\Bn 



g(rB £ r ) 



dz 



a) 2 $*V ± $ - 2flfl' <f>*d s <S> 
^ 2 $Vif + 2imfin' 



V 2 * 



(A.ll) 



The only new terms compared to equation [I?] are the 
ones containing the function g. One easily checks, by 
expanding g and £ r , that these terms are also hermitian. 
Therefore A.ll is also a variational form. We can draw 
the same conclusion as in the constant-_B case. 



A. 3. Dispersion Relation and Alfven flux 

Taking into account the gradient of B Q makes the analyt- 
ical derivation of the dispersion relation and Alfven flux 
much more complex and we will not attempt it here. 

If we make the additional assumption that Bq depends 
on r only weakly, i.e. that the current j° is weak, we can 
get rid of the influence of the terms containing g . This 
assumption is equivalent to requiring that the magnetic 
stress term is at most of the order of the centrifugal force, 
which should be the case in a realistic jet model including 
an equilibrium flow along the field lines. Using Vaoc as the 
Alfven velocity in the corona we obtain the condition: 
V (s) = [BoFJoc ~ r 2 n 2 

In this manner the conclusion presented in the body of 
this article for a constant vertical B° field can be extended 
to a weakly varying B®. 
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